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Parallel transport along circular orbits in orthogonally transitive stationary axisymmet- 
ric spacetimes is described explicitly relative to Lie transport in terms of the electric 
and magnetic parts of the induced connection. The influence of both the gravitoelec- 
tromagnetic fields associated with the zero angular momentum observers and of the 
Frenet-Serret parameters of these orbits as a function of their angular velocity is seen on 
the behavior of parallel transport through its representation as a parameter-dependent 
Lorcntz transformation between these two inner-product preserving transports which is 
generated by the induced connection. This extends the analysis of parallel transport in 
the equatorial plane of the Kerr spacetime to the entire spacetime outside the black 
hole horizon, and helps give an intuitive picture of how competing "central attraction 
forces" and centripetal accelerations contribute with gravitomagnetic effects to explain 
the behavior of the 4-acceleration of circular orbits in that spacetime. 
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1. Introduction 

Orthogonally transitive stationary axisymmetric spacetimes play a central role in 
our ideas about rotation and angular momentum in the context of general rela- 
tivity. Black hole spacetimes, the Godel spacetime, and even Minkowski spacetime 
expressed in coordinates adapted to this symmetry have shown how interesting and 
different rotation can be in the relativistic setting compared to our intuition from 
nonrelativistic mechanics and Newtonian gravitation. The constant speed circular 
orbits in these spacetimes are crucial in almost every coordinate system used to 
represent them, and their geometry reveals exactly how relativistic effects entangle 
matters associated with rotating observers. 

A family of test observers filling a region of spacetime and each individually 
rotating at a constant angular velocity (a circular observer family) is described by 
a "quasi-Killing" congruence, 1 namely a congruence associated with a vector field 
which is a stationary axially symmetric position-dependent linear combination of 
the two Killing vector fields describing the stationary axial symmetry of the space- 
time. For spacetimes which are asymptotically flat away from the symmetry axis, 
two privileged families of observers exist: 2 1) the static observers, or distantly nonro- 
tating observers, following the trajectories of the timelike Killing vector field which 
is regular in that limit but, in general, is characterized by a nonzero vorticity; 2) 
the zero-angular- momentum observers (ZAMOs), or locally nonrotating observers, 
which have zero vorticity and are orthogonal to the axisymmetry Killing vector 
which has closed circular orbits. The congruence of static observers determines the 
time coordinate lines while the family of orthogonal hypersurfaces of the ZAMOs 
determines the time coordinate hypersurfaces of the natural Boyer-Lindquist-likc 
coordinate system so often used in representing these spacetimes. 3 

The world lines of the ZAMOs and static observers respectively have zero values 
of the angular momentum "p^" associated with the axisymmetry Killing vector field 
and of the "angular velocity" <f> associated with the Killing vector field of the station- 
arity which is regular at spatial infinity away from the axis of symmetry, revealing 
two different aspects of the "nonrotating" property for orbits in the nonrelativistic 
theory which become distinct from each other in the relativistic theory when the 
spacetime itself is "rotating." In contrast to these two "extrinsic scalar quantities" 
for circular orbits which are defined in relation to these two Killing fields, one also 
has two intrinsic scalar quantities whose extreme values characterize another pair 
of "nonrotating" properties for a world line in the nonrelativistic theory, but which 
also become distinct in the relativistic theory in a rotating spacetime. These intrin- 
sic scalars are the Frenet-Serret curvature (magnitude of the 4-acceleration) and 
the magnitude of the Frenet-Serret angular velocity. 1 

At large distances the extremely accelerated observers 4-6 in some vague sense 
require the maximum outward acceleration to resist falling into a rotating central 
mass since they have the least centrifugal force needed to counteract the attraction, 
while the minimum intrinsic rotation observers (MIROs) 7-9 are rotating the least 
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in the sense that their Frenet-Serret frame, which is in turn rigidly related to any 
stationary axisymmetric frame for circular orbits, rotates the least with respect to 
local gyro-fixed axes among all circular orbits of different angular velocities at a 
given radius and angular location. Thus the usual nonrotating stationary circular 
orbits of nonrelativistic Newtonian theory (and of the static relativistic case) sepa- 
rate into 4 different kinds of orbits in the relativistic case of a rotating (stationary 
but not static) spacetime, each of which embodies a different aspect of the prop- 
erties we associate with nonrotating motion in nonrelativistic theory. The MIROs 
turn out to be important in describing the geometry of the acceleration and intrinsic 
rotation of the family of circular orbits, discussed in appendix A. 

The differential properties of a spacetime with symmetry, arising from the action 
of the symmetry group on the metric, arc bound together along a Killing trajec- 
tory by a continuous family of Lorentz transformations of the tangent space along 
it which relates parallel transport to Lie transport, since both transports preserve 
inner products. Given a frame of vector fields invariant under the symmetry group 
action, the induced connection matrix along a Killing trajectory, namely the value 
of the matrix-valued connection 1-form on the tangent vector to the curve, is an- 
tisymmetric when index-lowered with the metric reflecting the invariance of inner 
products, and generates this family of Lorentz transformations, modulo a constant 
linear transformation required to orthonormalizc the frame. Killing trajectories play 
a key role in interpreting the geometry of spacetimes which admit them and the 
geometrical properties of such curves associated with parallel and Lie transport are 
completely determined by this map, so it is of interest to understand its geometry. 

In particular, much has been written about stationary axisymmetric spacetimes 
and especially black hole spacetimes because of their importance in describing phys- 
ical systems and as models for understanding how rotation works in general rela- 
tivity. Circular orbits in these spacetimes are Killing trajectories which are cither 
helices or closed circles in the intrinsically flat world tube cylinders of the sym- 
metry orbits. These are the closest one can come to the Newtonian circular orbits 
in which our experience is grounded and have generated a great deal of interest 
over the past half century since the Godcl and Kerr solutions were found. While 
the usefulness of the literature on this topic is somewhat uneven, nevertheless its 
richness still offers interesting material to discuss, and the present problem seems 
worthy of exposition. 

The Serret-Frenet approach 1 ' 9 - 10 describes parallel transport along these Killing 
trajectory curves (when nonnull) representing constant speed circular orbits using 
a natural frame adapted to the local splitting of spacetime determined by their unit 
tangent vector U. On the other hand the language of gravitoelectromagnetism 11-13 
describes parallel transport in terms of the local splitting of spacetime associated 
with the 4-velocity u of either the ZAMOs or the static observers, both of whose 
world lines belong to the same larger family of all stationary circular orbits. The 
relationship between these descriptions involves the transformations of gravitoelec- 
tromagnetism, which show how the kinematical decomposition of the complete con- 
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nection along the observer family orbits adapted to their local splitting enters into 
the decomposition of the induced connection along a general orbit adapted to its 
own local splitting. In fact the whole electromagnetic analogy comes out of the 
decomposition of the covariant derivative of just the 4- velocity of a world line with 
respect to a family of observers. Here we extend the splitting to the entire covari- 
ant derivative along the world line for the case of transverse relative acceleration 14 
that occurs for circular orbits, showing how the gravitoelectric, gravitomagnetic 
and space curvature fields affect general parallel transport. 

A previous article, 15 following up work by Rothman et al 16 and by Maartens 
et al, 17 investigated parallel transport in stationary axisymmetric spacetimes along 
circular orbits which are confined to a plane containing oppositely-rotating circu- 
lar geodesies, like the equatorial plane of the Kerr spacetime or when additional 
cylindrical symmetry is present without translational motion in the additional sym- 
metry direction. This allows a study of clock effects, i.e., comparing periods of 
counter-revolving orbits defined in various ways, as well as circular holonomy. Par- 
allel transport of an initial vector along such an orbit leads to a variable rotation or 
boost in an invariant stationary axisymmetric 2-plane relative to the corresponding 
Lie dragged vector. The equatorial plane of the Taub-Nut spacetime 19 has more 
general behavior that is typical of the general case studied here, with both a boost 
and simultaneous rotation in a pair of mutually orthogonal invariant stationary ax- 
isymmetric 2-planes relative to the corresponding Lie dragged vector. The present 
article examines the geometry of parallel transport along (in general accelerated) 
nonnull stationary circular orbits in general orthogonally-transitive stationary ax- 
isymmetric spacetimes, appropriate to describe such orbits off the equatorial plane 
in the Kerr spacetime. Parallel transport then generically leads to a more general 
"Lorentz 4-screw" Lorentz transformation 20 generated by a 2-form whose interpre- 
tation in terms of electromagnetism helps put it into its canonical form correspond- 
ing to aligned electric and magnetic fields, a relationship not mentioned in current 
textbooks on relativity. 

The Kerr spacetime, originally studied very briefly in this context, 21 is used as 
an explicit example to illustrate these ideas. By constructing the explicit parallel 
transport transformation, the study of the Serret-Frenet properties of circular orbits 
begun in [9] is completed, and some some additional light is shed on the comparable 
strengths of the centripetal, gravitoelectric and gravitomagnetic effects along them. 

2. Induced connection on Killing trajectories 

The term "circular orbit" as used here refers to a Killing trajectory in a station- 
ary axisymmetric spacetime, i.e., a parametrized curve whose tangent vector is a 
constant linear combination of the two independent Killing vector fields generating 
the symmetry group. They may also be characterized as integral curves of Killing 
vector fields. General considerations as described in this section apply to any kind 
of symmetry group acting as isometries on spacetime. 
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Let e a (a = 0, 1, 2, 3) be a frame which is invariant under the symmetry group 
of a spacetime with metric g (signature — h ++), and let cu a be its dual frame. The 
connection components and connection 1-form matrix are 

V ea e / , = rT Q/3 e 7 , u} a p=T a lP u)i . (1) 

Given a parametrization of the Killing trajectory x a (X) — x a (c(X)) and the 
corresponding tangent vector c'(\) a — dx a (X)/dX, a family of Lorentz transforma- 
tions is generated by a fixed element of the Lie algebra of the Lorentz group which 
arises from the induced connection along the trajectory, namely the value of the 
connection 1-form matrix on the tangent vector of the parametrized curve, which is 
also Lie dragged along the curve (i.e., has a constant component matrix). If X is a 
tangent vector undergoing parallel transport along the curve, then its components 
in the invariant frame satisfy the constant coefficient system of linear first order 
differential equations 

» = T = ^+«^' < 2 > 

or 

dX a 



d , A a p X^ , A a Si = —w(c'(\)) a p = -PVdzVdA . (3) 

Since inner products are preserved by parallel transport, the induced connection 
matrix is antisymmetric when index-lowered to its fully covariant form 

0=— (e a - e/ 5) = -2A (a/3) . (4) 

Thus if one thinks of the constant matrix A = (A a p) as a mixed-tensor A = 
A a /3 e a ® oo 13 on the tangent space, this condition makes its fully covariant form 
A b = A a p uj a ® u 13 = \A a p uj a Alo 13 a 2-form and places the mixed tensor in the 
Lie algebra of the action of the Lorentz group of each tangent space along the 
curve. In an orthonormal invariant frame, the matrix then belongs to the matrix 
Lie algebra of the Lorentz group. Furthermore if one performs any change of frame 
e a — > e a i = (M^ 1 ) 13 a ep where dM^ a /d\ = holds along the Killing trajectory, 
then in the new matrix of induced connection, the coefficients are simply the new 
components of this mixed tensor A since there are no additional inhomogeneous 
terms in the transformation law for the induced connection along the curve. This is 
true for any new stationary axisymmetric frame along a circular orbit in the context 
of the class of spacetimes of interest. 

A vector which is Lie transported along such a curve has constant components 
with respect to the symmetry-invariant frame, so its covariant derivative is then 

Eigenvectors of the matrix A a p with eigenvalue zero therefore correspond to parallel 
transported fields along the curve which arc values of a single symmetry-invariant 
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vector field. When the tangent to the curve itself is such an eigenvector, the curve 
is a geodesic. 

The parallel transport system of equations has the exponential solution 

X a (\) = [e XA ] a XP(O) (6) 

representing a family of Lorentz transformations (modulo a constant linear trans- 
formation which orthonormalizes the frame). The nature of this family depends 
on the eigenvalue properties of the matrix generator A, which can only fall into 
one among four possible cases that can be classified by using the invariants of A 
under linear transformations. This task is performed by first introducing the de- 
composition of its index-lowered 2-form A b into electric and magnetic parts with 
respect to any future-pointing timelike unit vector u (representing the 4- velocity of 
an observer) 

A b = v 9 A £{u) + *WB(u) , (7) 

where *<- u >B(u) is the spatial dual in the local rest space orthogonal to u and £{u), 
B(u) are both 1-forms. While this decomposition is observer-dependent, under a 
change of observer these two fields undergo an electromagnetic boost under which 
the well-known quantities 

h=£(u)-B(u) = *A a pA^/A, I 2 = \\E(u)\\ 2 -\\B(u)\\ 2 = -A aP A^/2 (8) 

are Lorentz invariants related to the obvious matrix notation invariants through 
the relations det(A) = —l\ and Tr(A 2 ) = 2I 2 - It is convenient to refer to curves 
as electric-dominated or boost-dominated if I 2 > and magnetic-dominated or 
rotation-dominated if I 2 < 0, separated by the case I 2 = 0. 

In the terminology of Synge, 15,20 either A is nonsingular with Ji ^ (with 
four distinct nonvanishing eigenvalues) and a nonvanishing determinant (so A is 
also nonsingular in the ordinary matrix sense), or semi-singular with I\ = and 
I 2 7^ or singular with I\ — and I 2 — 0, the latter two cases having a vanishing 
determinant (so A is also singular in the ordinary matrix sense). The first case which 
is general corresponds to the generator of a simultaneous boost and rotation in a pair 
of mutually orthogonal 2-planes, the second case to the generator of a pure boost or 
rotation in a single 2-plane whose orthogonal 2-plane is invariant, and the last case 
to a null rotation when A is nonzero. The existence of a parallel transported vector 
(including geodesies, where the tangent vector is parallel transported) requires A 
to have a zero eigenvalue and so can only occur if A is at least semi-singular. 

For a nonnull Killing trajectory one can introduce the unit tangent vector U = 
c'(A)/|c'(A)-c'(A)| 1 / 2 , and use an arclength parametrization A = s, so U a = dx a /ds. 
One then has the decomposition 

A b = U b A£(U) +*<V)B(U) , (9) 

where the index-raised vector £([/)" = —a(U) = —DU/ds is the sign-reversed sec- 
ond (covariant) derivative of the parametrized curve. Its magnitude k = \ \a(U)\\ > 
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is the Frenet-Serret curvature, while £?([/)" is the Frenet-Serret angular velocity 
vector. To avoid complicating the discussion with additional signs, consider only 
timelike or spacelike 4- velocities U satisfying u-U < (on the future side of the local 
rest space LRS U of u within the whole tangent space), which can be parametrized 
by the relative velocity v with respect to the observer u 

U = f{u + ve) , 7= |l-t/ 2 | 1/2 , z/^±l (10) 

where e is a unit spacelike vector orthogonal to u, or in the corresponding null case, 
an unnormalized tangent 

v = ±1 : P± = u±e . (11) 

The limit [v, r )v 1 7] — > [±00, ±1,0] yields the Killing trajectories orthogonal to u. 
When U is timelike and future-pointing, then the electric and magnetic parts of 
the induced connection undergo a boost in the plane of u and U from (7) to (9) 
exactly like the electric and magnetic fields in electromagnetism. 11 The explicit 
construction of the exponential solution (6) of the parallel transport equations 
requires transforming A to a canonical form which allows the exponentiation to be 
performed explicitly, the canonical form being block-diagonal in the case of a non- 
null curve (closely related to the eigenvector diagonalization) , which is equivalent 
to Lorentz transforming a nonnull electromagnetic field to a frame in which the 
electric and magnetic fields are parallel. 

3. Circular orbits in a stationary axisymmetric spacetime 

Now specialize the discussion to stationary circular orbits in an orthogonally- 
transitive stationary axisymmetric spacetime. In Boyer-Lindquist-like (symmetry- 
adapted) coordinates (a; , x 1 , x 2 , x 3 ) — (t, r, 0, </>) the spacetime line element may 
be expressed in the form 

ds2 = d 4,0)+ ds M)' ( 12 ) 

where 

ds (t,4>) = 9tt dt 2 + 2gt0 dt Acj) + Acj) 2 , 

ds M) = 9rr dr 2 + g ge d6» 2 , (13) 

and the metric coefficients do not depend on the coordinates t and 4>. Let — 
g^ 2 di (i = r, 8, <p) be the orthonormalized spatial frame tangent to the constant 
time coordinate hypersurfaces. Orthogonal transitivity is just the condition that 
the directions along the orbits, namely the t-<f> 2-plane in the tangent space, are 
orthogonal to the r-6 2-plane in the tangent space. As will be seen shortly, it is 
natural to call these two subspaces the circular velocity plane and the circular 
acceleration plane respectively. 

The unit normal eg = n to the time coordinate hypersurfaces completes this 
spatial triad to a stationary axisymmetric orthonormal spacetime frame e&, and 
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represents the 4- velocity of the family of ZAMOs, a family of accelerated observers 
with 4-acceleration a(n). They are locally nonrotating in the sense that their vortic- 
ity uj(n) vanishes but they have nonzero expansion tensor 9{n) (minus the extrinsic 
curvature of the constant t hypersurfaces) whose nonzero components can be com- 
pletely described by the expansion vector #| = 0{n) a p due to the form (13) of 
the metric. 

Introducing the lapse function N = (— g") -1 / 2 and shift vector field with its 
only nonvanishing component = gt<j> / g<j><j> one has 

n=^[dt-N+dd, n b = -JVdt (14) 

and the metric itself is 

ds 2 = -N 2 dt 2 + g^ (<ty + N+db) 2 + dsf r g) . (15) 

The ZAMO kinematical quantities then have nonzero components only in the r-6 
2-plane of the tangent space 

a(n) = {a(nf, a(n) § ) = (d{\nN)/df, d(ln N)/d8) , 

H = \V9U (dN^/dr, dN*/86) , (16) 

where it is convenient to adopt a 2-vector notation indicated by an arrow for pairs 
of orthonormal frame components of vectors belonging to this subspace. Since the 
expansion scalar (trace of the expansion tensor) is therefore zero, the expansion 
tensor coincides with the shear tensor. In the static limit — > 0, the shear vector 
9^ vanishes. 

A general Killing trajectory is a stationary circular orbit with constant relative 
velocity along the Killing angular direction e^. The family of nonnull Killing vec- 
tor unit tangents U of uniformly rotating nonnull orbits at a given radius can be 
parametrized directly by the relative velocity v or by the rapidity a as follows 



{ 



U = 1 (±n + ve i ) (17) 

± cosha n 4- sinha , timelike: \v\ < 1, v = tanha , 
sinha n± cosha , spacelike: \i/\ > 1, v = ±cotha, 

where 7 = |1 — ^ 2 | -1 / 2 and in the future-pointing (+) timclikc case which will be 
of primary interest here: 7 = cosha. Note that for either case above 

dv_ _ ( 7~ 2 , U timelike, , , 

da \ =F7~ 2 , U spacelike. 
The angular velocity £ also parametrizes the half-family U ■ n < for ( e (—00, 00) 

u = r [dt + Cd<t>] , r = j/n, c = — + -— v . (19) 



It is convenient to refer to the cases \v\ < 1 and \v\ > 1 respectively as subluminal 
and superluminal relative velocities and extend the terminology of acceleration to 
the superluminal case. 
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These orbits are all open helices in spacetime except for the limiting case v~ x — > 
(or equivalently £ -1 — > 0) of closed (^-coordinate circles. For the null cases v = ±1, 
one must give up normalizing the tangent vector, letting instead P± = n ± e^. As 
long as ( ^ 0, one can use A = as a parameter along each such orbit, as done in 
[15] where loops in <f> were relevant to clock effects, or one can use the arclength 
s (proper time tjj) as a parameter in the nonnull (timelike) case, as is convenient 
for the Frenet-Serret approach. To translate the induced connection matrix defined 
on the circular orbit by Eq. (3) from one parametrization to the other for future- 
pointing timelike U with Q ^ 0, one needs the chain rule 

^ = ^2* Cr 5£ wherc CT^g-^j^-N^/N), (20) 

dru dru d<f> s d<f> s y ,y ' ' K 1 

so that 

^4(r(/-parametrization) = ( T yl(0-parametrization) . (21) 

The proper time (or arclength) parametrization will be used exclusively in this 
article in place of the ^-parametrization used previously. 15 

The family of future-pointing timelike 4-velocities is a branch of a hyperbola 
in the relative observer plane of n and e?, namely the t-(f> 2-plane in the tangent 
space spanned by the Killing vector fields, equivalently the subspace tangent to 
the cylindrical orbits of the symmetry group, called the circular velocity plane 
above. Another hyperbola in this plane describes the spacelike orbits, with tachy- 
onic 4-velocities, while the common asymptotes themselves correspond to the two 
oppositely rotating circular photons where no invariant preferred parametrization 
exists to normalize the tangent vector. 



4. The acceleration plane 

The corresponding curve of acceleration vectors for the family of future-pointing 
timclike circular orbits is instead a branch of a hyperbola 7-9 in the orthogonal tan- 
gent plane of the nonignorable coordinates r and 6>, called the circular acceleration 
plane above but hereafter called the acceleration plane for short 

a(U) = ^ = K(cos X ef + sin X e § ) a(U) = (a(U) f , a(Uf) , (22) 

where (k, x) are polar coordinates in the e r -eg plane in the tangent space, both 
constant for a given orbit. The Frenet-Serret angular velocity traces out a branch 
of a complementary hyperbola centered at the origin with the same asymptote 
orientations. The vertex of both hyperbolas occurs at the same relative velocity 
j/( vcrt -) which characterizes the MIRO family. This geometry of the acceleration plane 
is developed further in appendix A. The second branches of these two hyperbolas 
correspond to the spacelike orbits, for which it is convenient to extend the use of 
the word acceleration to refer to a(U), while the two asymptotes correspond to the 
acceleration vectors of the two oppositely rotating circular photon orbits. 
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Tabic 1. Geometric quantities for the Kerr spacetime, using the abbreviations A = 
r 2 -2mr + a 2 , A = (r 2 + a 2 ) 2 - a 2 s 2 A, S = r 2 + a 2 c 2 , (c, s) = (cos0, sin0). The 
orthonormal components of vectors in the r-9 plane are denoted by X = (X r , X 9 ). 

(JV, JV*) = ([AS/A] 1/2 , -2amr/A) 
(flrr, gee, gu) = (S/A, E, As 2 /E) 

3 („) = _|L- ([( r 2 _ a 2 c 2 )( r 2 + a 2 ) 2 _ 4rorW] A" 1 / 2 , -2ra 2 SC (r 2 + a 2 )) 
0. = (_ a s [ r 2 ( 3 r 2 + a 2 ) + a 2 c 2 (r 2 - a 2 )] , 2r a 3 s 2 c A 1 / 2 ) 

-- [(r 2 + a 2 )T? + 2mr a 2 s 2 {r 2 + a 2 + 



The acceleration a(J7) of the family of future-pointing timelike circular orbits 
in a stationary axisymmetric spacetime can be expressed in terms of the ZAMO 
relative observer decomposition in the following form (see equation (6.1) of [9], 
where the notation A is used for a(n)) 

a(U)=j 2 [fc (lic) v 2 + 26^ + a{n)] 

= fc(i ie ) sinh 2 a + 29^ sinh a cosh a + a(n) cosh 2 a 

1 - - 1 

= 2^ k ( lic ) + a(n))cosh2a + 9^ sinh2a+ -(-fc(i ie ) + a(n)) 

= - (a+ + a_ ) cosh 2a + - (a+ — a_ ) sinh 2a + a , (23) 

where fc(i; e ) are the components of the ZAMO Lie relative curvature vector fc(ii e ,c/,n.) 
(for more details, see [8,13]), and a± = fc( lic )+a(n)±2 9^ are the corotating/counter- 
rotating photon circular orbit accelerations (corresponding to a unit ZAMO energy 
affine parameter, ie., using P = n ± in place of U) which are aligned with the 
asymptotes of the acceleration hyperbola and ao = ^(a(n) — fc(H e )) is the center 
of the acceleration hyperbola. Fig. 2 of [9] shows some typical examples of these 
hyperbolas for the Kerr spacetime. 

Table 1 shows the geometric quantities needed to analyze the acceleration plane 
geometry for the Kerr spacetime. Discussion below will be limited to the case < 
a < m and to radii outside the black hole outer horizon at r( hor ) = m+(m 2 — a 2 ) 1 / 2 
where A = 0, with a > characterizing a black hole corotating in the positive <j) 
direction. In the equatorial plane 9 = ir/2, where parallel transport was discussed 
extensively in [15], the only nonvanishing components are radial 

rnA" 1 ' 2 [(r 2 |s 2 ) 2 -4a 2 mr] 
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ma(3r 2 + a 2 ) 

Si. : 



r 2 (r 3 + a 2 r + 2m a 2 ) 
_ A'/ 2 (r 3 -ma 2 ) 

(lie)r - r 2( r 3 + a 2 r + 2ma 2) ' I > 

The roots v+ > and i/_ < of the quadratic expression in v in square brackets 
in the first expression (23) for the purely radial acceleration in the equatorial plane 
define the ZAMO relative velocities of the corotating and counterrotating circular 
geodesies respectively (assuming a > 0) 

_ r 2 + a 2 =F 2 a (mr) 1 / 2 
V± ~ A 1 / 2 [a±r(r/m) 1 / 2 } ' 1 ' 

which are related to the quadratic form coefficients by 

a(n)f = fc (lic) . v + v_ , 6^. = -fc ( iic) f (y+ + v~)/2 • (26) 

In the context of spacetimc splitting techniques ( "gravitoelectromag- 
netism" 11 - 13 ), one can mimic the Newtonian situation by introducing a "relative" 
gravitational force as seen by the ZAMOs which is analogous to the Lorentz force of 
electromagnetism. Since n is a stationary vector field, its intrinsic derivative along 
U (see Eq. (5)) is 

F$ ) a = -^- = A° fi n°=£( n r (27) 

which lies in the acceleration plane and in turn may be expressed in terms of the 
ZAMO gravitoelectromagnetic vector fields by 

45 = 7 [9 - = 7 [9 + ^ x n H] , (28) 

where g = —a(n) are the components of the gravitoelectric vector field and H are 
the components of the effective gravitomagnetic vector field H = x„ 6^ due 
to the "differential rotation" of the ZAMOs arising from the shear of the ZAMO 
4-velocity field n. The latter field satisfies 0^ — —e^ x n H and \\9^\\ — \ \H\\, where 
x„ is the natural spatial cross-product in the local rest space of the ZAMOs. 

Similarly the covariant derivative of the stationary vector field e^, which defines 
the Fermi- Walker relative curvature vector, 22 leads to 

^ = 7 ^ fw) = -A^e? = -[*^B(n)r fj e/ = [B(n) x„ e $ ] a ■ (29) 
Since B(n) has no components along e^, this inverts to 

B(n)« = x„ [ 7 ^fc(f w )] or *™B(n) = Ja [ 7 ^%w)] b • (30) 

One can interpret this as the Fermi- Walker relative angular velocity vector asso- 
ciated with the changing direction of the trajectory in the local rest space of the 
family of ZAMOs along the world line (dividing out the gamma factor corresponds 
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to ZAMO proper time along the trajectory). Given these relationships, the accel- 
eration of the orbit can also be written in various ways 

a(U) = - 7 2 \g-2u0 $ - v 2 fc (lie) ] = - 7 + (7 vf fc (fw) 

= - 1 [£{nf +ue^x n B{nf] . (31) 

The Fermi- Walker relative curvature is related to the Lie relative curvature for 
a stationary circular orbit by Eq. (8.7) of [22] 

%w) = fc(iic) + • (32) 

Using this to rc-cxprcss B then leads to the final pair of equations 

£{nf = 1 [g + ue^ > x n H] = 7 [g - v6£[ , 

B{nf = 7^ x„ fc (lie) + H]=je^x n [v k (lic) + 6$] . (33) 

For the ZAMO orbits {v = 0, 7 = 1), the electric part of A is just the gravitoelectric 
vector field and the magnetic part of A is just the effective gravitomagnetic vector 
field. 

The relations (33) almost have the form of an electromagnetic boost except for 
the two distinct accelerations, but introducing the sum and difference fixes this 
leaving an extra correction. Letting g± = (g ± fcn ie ))/2, one finds 

(34) 

where 

, ,. . / 7 7^e;X„\ / cosha sinhaej-x^X ,„_. 
M(a) = ' = . 1 , 35) 

\--fiye^x n 7 J \-sinhae^x„ cosha J 

is the boost operator for a pair of electric and magnetic fields in the r-6 2-planc 
LRS n n LRSu orthogonal to the direction of relative motion in the local rest space 
subspace common to both n and U, where the projection P(n, U)^ 1 of the general 
formula (4.14) of [11] acts as the identity transformation. This is true for world lines 
characterized by purely transverse relative acceleration, 14 as are circular orbits seen 
by circularly orbiting observer families. The rapidity parametrization of the boost 
matrix here only holds when U is future-pointing and timelike; the remaining cases 
of Eq. (17) must be considered separately. 



5. Parallel transport along circular orbits and the Frenet-Serret 
description 

The index-lowered matrix governing parallel transport along an arclcngth (proper 
time when timelike) parametrized circular orbit with 4-velocity U is 

A a/3 = -{TQ [flV>[a,/3] + C 1 9t[a.,l3]] = -^ M ff M [a,/3] (36) 

where the factor T( — 7C/A = d^/ds takes into account the arclength / proper 
time parametrization of the trajectory and would instead have the value 1 if the 
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curve were parametrized by the angle 4> as in previous articles. 15, 19 The 2-form A b 
can be expressed in terms of its electric and magnetic parts either with respect to 
U or with respect to the ZAMOs 

A b = U b A £{U) + *(U)B(U) = n b A £{n) + *^B(n) . (37) 

Because of the orthogonally-transitive stationary axisymmetric form (13) of the 
metric components, the matrix (A a p) can only have nonvanishing components for 
index pairs (t, r), (t, 6), (4>, r), (<fi, 9) in either order. This means that its correspond- 
ing electric and magnetic vector parts can only have nonzero components in the 
r-6 plane, namely the acceleration plane where all the acceleration terms lie. The 
transformation between the electric and magnetic quantities relative to U and those 
corresponding to n defines the electromagnetic boost already introduced above 

(£(U)*\ (-y[£{n)* + ve t x n B(n)t\\_ (£(n)*\ 

Thus this represents a transformation of electric and magnetic fields in the acceler- 
ation plane as well. Composing this with (34) leads to the fact that the U fields only 
involve the relative velocity through an inhomogeneous boost by twice the rapidity 
of the first boost 

/■ £ <^-M (2Q > (■'-) + (39) 



\B(uy J v ; \H J \ 

recalling that g± = (g±k(u e A /2. The extra term for the electric part g + = —So is the 
sign-reversed center of the acceleration hyperbola, representing the shift of origin 
needed to follow the double boost which takes the (modified) ZAMO gravitoelectric 
and gravitomagnetic vectors to the electric and magnetic parts of A. 

The decomposition of A with respect to U can be interpreted by expressing it 
in terms of the associated spacetime Frenet-Serret frame {e = U,e\,e2,ez\ and 
the curvature and torsions k, t\ and t 2 1,9 

De De l De 2 De 3 

- — = net , - — = k e + ti e 2 , - — = -T\ e x + r 2 e 3 , - — = -r 2 e 2 . (40) 

OTu dTjJ CtTu OTu 

From the identification 

De 



dru 

one can read off the identifications 



= -AP a ep, (41) 



A^ = —k U A ei + T\ e\ A e 2 + t 2 e 2 A e 3 , 
£(U)* = -a(U) = -k ei , B(U)* = lu (fs) = n e 3 + r 2 e x , (42) 

assuming that {ei, e 2 , e 3 } is a right-handed triad (so that e 3 = e\ e 2 ). Both 
£(U)$ and B(U)$ belong to the acceleration plane span({ei, e 3 }) = span({ef , e§}), 
while {eo, e 2 } span the velocity plane. By continuity in the rapidity a, the Frenet- 
Serret frame may be easily extended to geodesic orbits where k = 0. 
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The Frenet-Serret scalars completely describe parallel transport along a circular 
orbit with respect to a frame determined by the curvature properties of the orbit 
itself. Because the unit tangent of the constant speed circular orbit is stationary, so 
too is its entire Frenet-Serret frame, which on a given trajectory is related to the 
natural symmetry adapted frame associated with the Boyer-Linquistlike coordinates 
by a constant linear transformation. 

The Synge classification of the matrix A is easily translated into conditions on 
the intrinsic Frenet-Serret quantities 

h = £(U) ■ B(U) = -a(U) ■ uj {fs) =-kt 2 , 

h = \\£{U)\\ 2 \\B(U)\\ 2 = « 2 ~ INfs)!! 2 = « 2 - (r! + r|) . (43) 

Thus A is nonsingular (I\ ^ 0) in both the Synge sense and the ordinary linear 
algebra matrix sense if the orbit U is nongcodcsic (k ^ 0) and has nonzero second 
torsion. A is semi-singular {I\ = but I 2 ^ 0) if either U is geodesic or the second 
torsion is zero, while it is singular (ii = = I2) if in addition Frenet-Serret angular 
velocity has the same magnitude as the acceleration. 

For orthogonally-transitive stationary axially symmetric spacetimes, the two 
torsions are simply related to the polar form of the acceleration vector in the ac- 
celeration plane under certain assumptions about the choice of the directions of 
the Frenet-Serret frame vectors 9 as described in detail in Appendix B. For the 
future-pointing timclike case, one has 

/TT-, , • \ 1 dn 1 dx .... 

a(U) = K{cosxe f +smxe 9 ) , n = - - — , r 2 = - - k — , (44) 

which can be converted back to ^-derivatives using Eq. (18). Vanishing first torsion 
corresponds to a critical point of k as a function of the rapidity a (including the case 
of the extremely accelerated orbits) , while vanishing second torsion corresponds to 
a critical point of the polar angle X- 

Finally, note that the corotating and counterrotating circular photon orbit tan- 
gent vectors P± — n±e^ are parallel transported along U when D(n±e^)/dTij = 0. 
This translates into the statement 

£{n) = Dn/rjj = ±De^/dT V = ±B(n) x„ , (45) 

which implies that £{n) and B{n) are orthogonal and equal in magnitude, so both 
invariants I\ = = I2 vanish corresponding to the Synge singular case. 

Fig. 1. The Frenet-Serret scalars plotted versus relative velocity v at r/m = 4, 3, 2 (top to 
bottom) and 8 = tt/2, it/4 (left to right) for future-pointing timelike circular orbits in an a/m = 
1/2 Kerr spacetime: k (solid), t\ (short-dashed), T2 (long-dashed, absent for 8 = tt/2 where T2 = 0) 
and ||w(fs)|| (dotted, just |ti| for 8 = tt/2). The radii chosen for these three figures belong to the 
equatorial plane radial intervals where a) two oppositely-rotating timclike geodesies exist, b) one 
corotating timelike geodesic exists, and c) no timelike geodesies exist, which for 8 = n/4 roughly 
correlate respectively with the regions A, B, C in Fig. 1 of [9] characterized by the same number 
of timclike orbits with zero radial acceleration at a given (r, 8) location. 
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a) r = 4m > r goo - , 9 = n/2 6 = 7r/4 
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To make these abstractions more concrete in the Kerr spacetime, the Frenet- 
Serret frame has been defined so that its limit as r — > oo has the slow-motion far- 
field Schwarzschild properties of circular orbits: the acceleration is approximately 
radially outward while the boosted direction of increasing (j) is approximately along 
the ^-coordinate direction, leading to 

ei ~ e f , e 2 ~ , e 3 = eiX[,e 2 ~ -e§ . (46) 

The leading contribution to the first torsion is the special relativistic angular ve- 
locity T\ ~ vjr of the rotating radial unit vector along a circular orbit of ZAMO 
relative velocity v. For a static spacetime this remains an odd function of v, but 
nonzero rotation introduces asymmetry, shifting the zero of t\ away from v = 0. 

Fig. 1 shows the typical behavior of the Frenet=Serret scalars plotted versus 
relative velocity both on the equatorial plane = ir/2 and off (9 — 7r/4) for an 
a/m =1/2 Kerr spacetime at radii approaching the black hole for the three radial 
intervals where two, one and then no timelike circular geodesies exist in the equato- 
rial plane. Figs, la) and 16) show a kink in the curvature k in the equatorial plane 
at its zeros (geodesies), with a resultant sign change jump in t\. By allowing the 
Frenet-Serret vector e\ to remain equal to in the equatorial plane, while allowing 
k to assume negative values there as done in Figs. 5 and 9 of [13], both graphs 
become diffentiable at the location of the geodesies where the Frenet-Serret frame 
is only defined by continuity with the nearby accelerated orbits. This sign change 
for the torsion makes its new graph coincide with its absolute value | |w(fs) 1 1 (dotted 
curve) in the current plot. 

In Fig. la) as one moves up off the equatorial plane at 9 — ir/2, the three 
critical points of the curvature k, all local extrema representing the two geodesies 
and the extremely accelerated orbit on that plane, all evolve, eventually resulting 
in only one critical point below a certain angle 9 representing the single zero of the 
first torsion. Off the equatorial plane where it is identically zero, t 2 is a monotonic 
function of v which has exactly one zero at the extreme value of the polar angle x 
parametrizing the acceleration. 

The fact that for sufficiently large positive v 1 t 2 is negative, and positive for 
sufficiently negative v shows that \ is respectively increasing and decreasing in 
those limits, so the side of the hyperbola closer to the "horizontal" dashed line 
in Fig. 2a) of [9] is the sufficiently positive v side of the hyperbola, while the 
sufficiently negative side is pushed farther away from the horizontal, both of which 
continue rotating away from that horizontal direction as one approaches the horizon 
in Figs. 26) and 2c). The true orientation of these hyperbolas is more apparent 
in the limiting Schwarzschild Fig. 5 of that article, where the chief affect of the 
gravitomagnctic field in the Kerr spacetime is to split the half-ray of gravitoelectric 
acceleration vectors into the two halves of the hyperbola branch, with the corotating 
half closer to the horizontal sufficiently far away from the horizon; the hyperbola 
rotates up and away from the vertical direction and outwards towards the radial 
direction as the horizon is approached. This is simply because the curvature vector 



0.15 



0.05 



-0.05 



Geometric transport along circular orbits in stationary axisymmetric spacetimes 17 

7 2 : 6> = 7r/2 Zi: 6> = 7r/4 




0.5 



0.25 



-0.25 



-0.5 





/ 87 / / //|| 
(^^^^7=2.87 







-1 



-0.5 



0.5 



1 



Fig. 2. The invariant I2 plotted versus relative velocity at 9 = ir/2 where Ii = for ana/m = 1/2 
Kerr spacetime at selected radii (the situation at 8 = 7r/4 is not much different qualitatively) and 
the invariant I\ at = n/4. 



fc( lic ) is well-behaved while the ZAMO acceleration a(n) (A in that figure) blows 
up as the horizon approaches and so dominates their linear combinations, and the 
latter acceleration field is approximately in the outward radial direction 9 to keep 
a test particle from falling towards the "center of attraction" when "at rest with 
respect to the ZAMOs." This is the reason the electric part of the connection must 
dominate the magnetic part as one approaches the horizon in the Kerr spacetime: 
modulo common gamma factors the centripetal effects associated with the dominant 
term in the magnetic part are bounded relative to the growing ZAMO acceleration 
which dominates the electric part (see Eq. (33)). 

Fig. 2 shows instead the two invariants at selected radii on and off the equato- 
rial plane. The graphs of I2 in the equatorial plane show how the boost zone where 
it is positive expands as one approaches the horizon, first going superluminal for 
counterrotating orbits and then later for corotating orbits. 72 and hence I\ is iden- 
tically zero in the equatorial plane, but off the equatorial plane, I\ is an increasing 
function which grows faster as one approaches the horizon, always passing through 
zero at a very small negative relative velocity not visibly different from zero at this 
scale where T2 = 0. 



6. The general solution of the transport equation along 
future-pointing timelike circular orbits 

As shown in section 2, the solution of the parallel transport equation is (6), where 
A is given by (36) and A = r is the proper time parametrization of the orbit when 
U is timelike. A is a tensor field only defined along the orbit and the components of 
A can be referred to any observer-adapted frame once the observer family itself is 
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chosen all along the orbit. Splitting A into its electric part associated with a boost 
in an observer-adapted frame and a magnetic part associated with a rotation 

A b = u b A £{u) + *™B{u) = A B (u) b + A R (u) b , (47) 

can be done with respect to any stationary unit timelike future-pointing vector field 
u representing the 4- velocity of a stationary observer family defined along the orbit, 
but the separation is observer-dependent. 

In the general Synge nonsingular case, A generates a simultaneous boost and 
spatial rotation acting in orthogonal 2-planes, but the frame is not adapted to 
these 2-planes so the boost and rotation are mixed together in their effects on the 
frame vectors. The two parts Ab(u) and Ar(u) do not in general commute, so 
e Ar ^ e A B (u)r e A K (u)r_ However, one can always find one observer for which these 
two parts do commute, which in turn corresponds to parallel electric and magnetic 
fields £{u) and B(u). This has been explained in detail in [3] (exercise 20.6, p. 480) 
for a nonnull electromagnetic field where a new observer with relative velocity along 
the direction of the Poynting vector can be chosen to align the electric and magnetic 
fields. We use that same procedure here. 

For simplicity assume U is a future-pointing timelike unit vector. Start with 
A b — U b A £{U) + *<-U)B(U) as seen by U itself. Pick a new circularly rotating 
observer with 4-velocity X(U) traveling in in the direction of £(U) Xjj B(U) with 
signed relative velocity tanhd along the <fi direction within the local rest space of 
U chosen so that 

As long as the right hand side does not equal 1 (the Synge singular case for A), 
this has a real solution. With a convenient choice of sign, direct evaluation using 
Eq. (42) gives 

2 K T1 

tanh 2 a = -z ,, = T . (49) 

« 2 + IKfs)II 2 

while identities then lead to 

1 - Vl - T 2 t 1 

tanh a = — = t , sinh a = , , cosh a = ■ . (50) 

Note that the sign of a is the same as the sign of t\ as long as k > as 
assumed here, while a itself reduces to in the limit n — > or n — > and the 
Frenet-Serret frame is itself adapted to the parallel transport boost and rotation. 
The first case of vanishing first torsion corresponds to the extremely accelerated 
orbits; when in addition the second torsion goes to zero as it does identically in the 
equatorial plane of the Kerr spacetime, this becomes a pure boost. The second case 
of vanishing curvature corresponds to the geodesies in the equatorial plane, which 
instead leads to a pure rotation in the local rest space of the geodesic. When a^O, 
its sign indicates whether one must slow down (— ) or speed up (+) in the positive 
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Fig. 3. The boost velocity tanhct from U to the parallel transport adapted frame plotted versus 
v, shown for r = 2, 2.5, 3, 3.4, 4 (from solid to increasingly longer dashed lines) in the a/m = 1/2 
Kerr spacetime. 



(^-direction relative to U to achieve the canonical form for the connection matrix 

Ay. 

Fig. 3 illustrates the behavior of the relative velocity tanha. Far enough from 
the horizon where there are two geodesies in the equatorial plane, a is zero at the 
two geodesies and at the extremely accelerated orbit in between them, but as one 
moves in, first the extremely accelerated orbit and the counterrotating geodesic 
go null together, leaving one zero, and then finally the corotating geodesic goes 
null (not yet reached in this sequence). Similarly the values ±1 occur when the 
Synge singular orbits are reached at the endpoints of the boost zone, which grows 
outwards from v = as one moves towards the horizon, exiting first at v = — 1 
from the timelike interval. Off the equatorial plane, there is only one zero at the 
extremely accelerated orbit. For velocities faster than this, one must slow down, 
and for velocities less than this, speed up, to boost U by this relative velocity. 

Boosting the Frenet-Serret frame vectors {U, e 2 } to {E n , E 2 } by the relative 
observer boost from U to X(U) and defining the second pair of frame vectors E\ 
and E 3 as proportional to their respective covariant derivatives gives the form 

Eq = X{U) = cosh a U + sinha e 2 , 

E 2 = X(U) = sinha U + cosh a e 2 , 

Ei = (Tg 1 [(k cosh a — ti sinhd) ei + r 2 sinhd 63] , 

E 3 = a^ 1 {(k sinha — n coshd) ei + r 2 coshd 63] , (51) 

(reducing to the identity transformation when a = 0) which leads to 

A b = -a B E b A E\ + c7r E\ A E\ . (52) 
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Fig. 4. <7g versus v for the a/ra = 1/2 Kerr spacetime. 
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Requiring orthogonality of E\ and E 3 leads to the condition (49), while the invari- 
ants (recall the relations (43)) 



-ob or 



-KT 2 , h 



4 



(53) 



may then be used to determine the nonnegative boost parameter <tb > and the 
rotation parameter or in terms of the curvature and torsions 



V = [(h 2 + h 2 ) 1/2 + h]/2 , ctr 2 = [{h 2 + h 2 ) 1 ' 2 h}/2 , 



(54) 



equivalent to normalizing E\ and E^. Note that E2 — X(U) is just the unit vector 
in the direction of increasing <f) in the local rest space of X(U) and the sign of <jr 
must agree with the sign of r 2 . 

In the Synge semi-singular case I\ = —KT2 = but I2 7^ 0, then parallel 
transport reduces to a pure boost (I 2 > — > a R = 0, a B = l^l 1 ^ 2 ) under which 
X(U) is invariant or a pure rotation (I2 < — > <t_b = 0, cjr = Ih] 1 ^ 2 ) under 
which X(U) is invariant. This is true in the equatorial plane of the Kerr spacetime 
where r 2 = identically and this invariant unit vector X(U) or X(U) respectively 
defines the equatorial plane map U — > Z(U) which picks out the direction of the 
covariant constant vector in the velocity plane along the orbit. 15 For geodesies in 
the equatorial plane, this is a pure rotation in the local rest space of U as remarked 
above. Figs. 4 and 5 show the typical behavior of <jb and /jr. 

To express the boosted frame in terms of the original symmetry adapted ZAMO 
frame, one can follow the same approach, starting with a boost 



Eq = X{U) = cosh an + sinhde^ 
and using its inverse to re-express 



E 2 = X(U) = sinha n + coshd e? (55) 



A k =n k Af(n)+^AX: , /C = 7 z/% w) 



(56) 
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Fig. 5. or versus v for the a/m = 1/2 Kerr spacetime. 

in the form 

A b = X(U) b A £(X(U)) + X(U) V A K{X(U)) , (57) 

where 

£ (X(U)) = cosha£(n) — sinhd/C , JC(X(U)) = — sinhd£ (n) + coshd/C . (58) 
Then the orthogonality condition (alignment of electric and magnetic vectors) 

= £(X(U)) ■ K(X{U)) = X(U) ■ [£(X(U)) x x{u) B(X(U))} (59) 
yields the corresponding rapidity condition 

£{n)-hZ 



tanh 2a 



• [£{n) x n B(n)} 



(60) 



,|£(n)|P + ||/C|| 2 ||f(n)|| 2 + ||B(n)|| 2 ' 

from which one can determine the ZAMO relative velocity vx = tanhd = [1 — (1 — 
tanh 2 2a) 1 / 2 ]/ tanh 2d of X(U) in terms of a(n), fc(n e ), and 9^ using (33). 

Now that the boost velocity is determined with respect to either the ZAMOs 
or U itself, consider the component matrix (A a p) in the new boosted frame E a , 
which has the canonical block-diagonal form 



(A a p) = (A B a p ) + (A R %) = 



<y b 
cr B 



2 



/0 2 
2 



2 

OR 

-<jr o 



(61) 



V o 2 o 2 y v L-^r o j / 

where A B has eigenvalues A = 0, 0, ±ctb and A R has eigenvalues A = 0, 0, ±it7R. 
This translates the action of parallel transport along the orbit for a generic vector 
X into two commuting actions of a family of boosts and rotations when expressed 
in this frame 



X(r)=AX(r) , 



X{t) = e A * T e ART X(0) , 



(62) 
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so that the frame components of a parallel transported vector X(t) in this new 
frame are only coupled within two mutually orthogonal subspaces 



X\r) 
X 3 (t) 



{ cosh erg t sinho'B''"\ /X°(0) 
y sinh (7b t cosh (7b t J \ V X 1 (0) 



cos ctr t sin ctr, t \ / X 2 (0) 
- sin (TrT cos or, t) \X 3 (0) 



(63) 



(64) 



and evolution of the initial data in the general case is governed respectively by a 
boost or a rotation in these subspaces. In the general nonsingular case, both must 
be present. 

When A is singular in the Synge sense, corresponding to vanishing values of the 
quantities cb and <7r defined by Eqs. (54) the above derivation and the represen- 
tation (52) do not hold and must be redone. Circular orbits for which A is singular 
but nonzero cannot be geodesies but must have nonzero k = ±ti and t 2 = 0. Then 
the mixed component matrix of A with respect to {u } e\, e 2} 63} is 



A b = K e 1 A[U±e 2 ], (A a p) 



/Ok 0\ 

K ±K 

±k 

\0 0/ 



(65) 



which generates a null rotation 



= 1 + At + ^-A 2 t 2 



(66) 



in the plane orthogonal to e$. 

In the Kerr spacetime, this Synge singular case only occurs in the equatorial 
plane as the transition between the parallel transport boost and rotation velocity 
zones for each circular orbit. 15 That A b is a decomposable 2-form for the Synge 
singular case is true for the entire Synge semi-singular case which characterizes 
the Kerr equatorial plane where the vanishing of t-i makes <tb Ei and (Jr E 3 both 
proportional to e\ (so only one term survives in Eq. (52), proportional to ei A Eq or 
ei A E 2 ), which then factors out of A^ as one also sees directly from setting t 2 = 
in Eq. (42) 



A* 



zU A e 



■i + ti ei A e 2 = ei A (k U + t\ e 2 ) . 



(67) 



When the vector in parentheses is nonnull, it is proportional to Eq = X(U) if 
timclikc (Ii > 0) corresponding to a parallel transport boost or E 2 = X(U) if 
spacelike (I 2 < 0) corresponding to a parallel transport rotation, and taking the 
paramctrization relationship (21) into account, while choosing instead e\ = ef 
(allowing k to change sign), one can translate this back to the map Z determining 
the invariant parallel transport direction using (A. 7) of [15]. 
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7. Boost and rotation domination zones 

The Synge class of the matrix A associated with the invariants I\ and I 2 is a 
function of the two nonignorablc coordinates r and 6. Either I2 ^ (nonsingular 
case) or I\ = but I 2 ^0 (semi-singular case) or I\ = = I 2 (singular case), each 
case corresponding to a different class of Lorentz transformations governing parallel 
transport. 

The sign of I2 indicates the relative magnitudes of the boost and rotation rates 
along a circular orbit in the nonsingular case (boost-dominated if I2 > 0, rotation- 
dominated if I2 < 0) and whether a pure boost (I2 > 0) or pure rotation (I 2 < 0) 
occurs in the semi-singular case. Such boost-dominated and rotation-dominated 
zones occur in each family of circular orbits at a given radius and angular location 
as a function of the relative velocity, but as one approaches the vicinity of the 
horizon in a Kerr spacetime, for timclike circular orbits the rotation-dominated 
zones disappear as the attraction to the black hole (acceleration) overcomes any 
relative centripetal acceleration effects (Frenet-Serret angular velocity) due to the 
limits on the relative velocity \i/\ < 1. 

For a given radial and angular position of a circular orbit, the condition I 2 = 
determines two relative velocities which delimit an interval V{pt-) < < V{pt+) 
around the ZAMO angular velocity v — within which parallel transport boost- 
dominance occurs (I 2 > 0) 

|*| > IM - ||£(n)|| > - > 7 ^ll%w)|| • (68) 

and the Fermi- Walker relative gravitational force on the circular orbit exceeds the 
Fermi- Walker relative angular velocity. Using (33) one can express this inequality 
as a quadratic inequality in v 

h = (| |^| | 2 - P (li c)||V 2 -26$- (g + k (]iB) )u+ \\g\\ 2 - ||^|| 2 

= -Av 2 - 2Bv + C > (69) 

whose corresponding equality has the solutions i^(pt-) < ^(pt+) given by the min- 
imum and maximum values of the usual quadratic formula for the roots 

-B±[B 2 + AC)}^ 
f(PT±) = ^ • (70) 

This leads to the static limit 

n llffll Vto + fcqic)^ , 7n 

01 -> : i>(pt±) -» jr, rr ±1 - | , rr- , (71) 

P(lio)ll \ |M| P(lic)ll J 

with B — 6^ ■ (g + fc(i; ) ) > for a corotating (a > 0) Kerr spacetime 9 and A = 
ll^(iic)l| 2 ~ ll^^ll 2 > holds everywhere except very close to the horizon. 

When the relative velocity moves out of the interval [^(pt-), ^(pt+)L tne circular 
orbit has sufficient speed for the relative Fermi- Walker angular velocity to exceed 
the Fermi- Walker spatial gravitational force. Inside this interval the Fermi- Walker 
spatial gravitational force dominates the relative angular velocity. At these endpoint 
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values, the circular orbit parallel transports one of the two circular null tangent 
vectors, as shown by Eq. (45); in the equatorial plane, A is then singular generating 
a null rotation as described in [15]. 

The condition ^(pt±) = ±1, after clearing fractions, isolating the radical and 
squaring both sides of the resulting equation, becomes A(A±2 B — C) = 0. Setting 
the second factor to zero yields 

(a(n) - fc (lic) ) • (a(n) + fc (lie) ± 2 9$) = <-» a ■ a± = (72) 

which shows that these limiting null orbits occur when the position vector ao = 
—g+ of the center of the acceleration hyperbola is perpendicular to one of the two 
asymptote directions a± corresponding to null orbits, or if a*o vanishes, or if one 
of those null circular orbit acceleration vectors vanishes. The latter corresponds to 
corotating and counterrotating null circular geodesies, which occur in the equatorial 
plane of the Kerr spacetime. 

The limits ^(pt±) — * ±°o correspond to S7pt±) = V^fPTi) — * ^ ± wmc h de- 
scribe the spacelike closed ^-coordinate circles. These occur where the leading co- 
efficient —A = H^H 2 — l|fc(iio)l| 2 g° es to zero in the above quadratic relationship 
(69) for the root where no cancellation occurs in the numerator in this limit. But 
for null relative velocities v = ±1, one has fy fw ) = fc( lic ) + (±1) 6^, so this will occur 
at locations (r, 0) where the Fermi- Walker curvature for one of the two null orbits 
vanishes, corresponding to a "Fermi- Walker relatively straight" curve, where the 
Fermi- Walker relative centripetal acceleration goes to zero. 

For a corotating Kerr spacetime, as one moves closer to the black hole in the 
equatorial plane, 15,23 the interval [k7pT-)i ^(pt+)] expands until z^(pt-) = — 1 at the 
counterrotating null circular geodesic so that counterrotating rotation-dominated 
timelike orbits no longer exist and then continues until finally ^(pt+) = 1 at the 
corotating null circular geodesic so that not even corotating rotation-dominated 
timclike orbits exist beyond that. 23 In the equatorial plane these occur exactly at the 
counterrotating (77) and corotating (ru) null geodesic orbits, while ^(pt-) - * — 00 
at the geodesic meeting point observer horizon (rm) and i^(pt+) 00 at the 
black hole horizon (ZAMO horizon: r^), corresponding to transport around closed 
<p coordinate loops. 23 Fig. 6 shows the boost-dominated zone (or just boost zone 
since there is no simultaneous rotation) shaded in grey, with the remaining rotation- 
dominated zone (or just rotation zone) in white. This explains the domains of <jb 
(inside the boost zone) and or (outside the boost zone) respectively in Figs. 4a) 
and 5a), while off the equatorial plane both quantities are defined for all velocities 
except \v\ — > 1 where they both go infinite due to a gamma factor in their definition. 

Fig. 7 extends these boundary regions off the equatorial plane to show the hori- 
zon, ergosphere boundary and the curves in an r-6 section of the exterior Kerr 
spacetime where z^pt±) take the values ±1 and have the limits ±00. The first two 
conditions define two closed surfaces of revolution surrounding the horizon and 
ergosphere boundary, the II surface (corotating) inside the I surface (counterro- 
tating), which both meet the horizon and ergosphere at the symmetry axis 8 = 
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Fig. 6. The boost- (grey shading) and rotation-dominated (white) relative velocity zones for the 
equatorial plane in the a/ra = 1/2 Kerr spacetime for — 00 < v < 00 using the reciprocal velocity 
v = 1/u to plot the spacelike relative velocities > 1, with the top and bottom borders identified 
to describe the same closed rf> circles being traced out in opposite directions. The velocity curves 
labeled by Z± = i^(pT±) separate these two zones. Vertical lines indicate the radii at which the 
countcrrotating (rj) and corotating (ru) geodesies go null, the radius rjn where the geodesic 
meeting point orbits go null and the horizon where the null meeting point orbits (ZAMO 
world lines outside the horizon, see [18]) go null. These correspond respectively to the radii where 
U (PT— ) = ~1 COj "(PT+) = 1 (-f-Oi "(pt-) ~* ~°° (HI) an d ^(PT-H) ~~ * 00 • Also shown are the 
relative velocity profiles of the geodesies (geo±), geodesic meeting point curves (gmp), and the 
critical points of the circular orbit acceleration where n = (crit±). The subluminal critical point 
curve (crit— ) corresponds to the extremely accelerated observers. The dashed curve between the 
radii rj and r/j shows instead the critical points of t\ = |o)(fs)I- 



but extend farther out from them and each other as one approaches the equatorial 
plane, which they intersect in circles at the respective radii where the corotating null 
geodesies and counterrotating null circular geodesies exist. The limit z^(pt-) — * —00 
(approaching from the side away from the origin, but 00 when approaching from 
inside) occurs along the smaller curve III in Fig. 7 near the equatorial plane up to 
about 9 = 71° where it intersects the horizon, where in turn the limit f(px+) — * 00 
occurs. 

In the static Schwarzschild limit, the surfaces corresponding to the curves I and 
II in Fig. 7 come together at the radius r — 3m in the equatorial plane, while the 
horizon, ergosphere boundary and the infinite limit surfaces collapse to the sphere 
r = 2 m (see Fig. 2 of 15). Inside the single surface / = II, no timelikc timelikc 
circular orbit exists with rotation-domination. In the general Kerr spacetime, in- 
side of the outermost surface /, no timelike counterrotating circular orbit exists 
with parallel transport rotation-domination, while inside the next surface II, no 
timelikc corotating or countcrrotating timelike circular orbit exists with rotation- 
domination. Inside the next surface III, no counterrotating circular orbit of any 
causality type exists with rotation-domination, while at the horizon all circular 
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Fig. 7. A polar plot in the r-0 plane treated as though they were polar coordinates in a flat 
2-space with horizontal axis p = (r/m)s'm9 and vertical axis z = (r/m) cos#, for an a/m = 1/2 
Kerr spacetime, showing the event horizon (h, thick curve), where i/( PT+ j — > oo, the ergosphcrc 
boundary (e, dashed curve), the curve / where i/( PT _j = —1 and the curve 77 where i^(pt+) = 1 , 
and finally the curve III intersecting the horizon at about 19° above the equatorial plane, along 
which ^(pt-) — * — oo from the outside, but oo from the inside (see Fig. 6). 

orbits of any kind are boost-dominated as the rotation-domination velocity zone 
shrinks to the null set at the closed (^-coordinate circle being traversed in the coro- 
tating direction of increasing <j>, exactly as in the equatorial plane in Fig. 6. 

One can imagine how this latter figure changes as one decreases 9 from 7r/2 in 
Fig. 6. Fig. 8a) shows how the equatorial plane velocity curves deform and join to 
produce the vanishing t\ velocity curves as the boost zone squeezes down towards 
v = and its marking radii move in towards the horizon as 9 decreases. This 
continues in b), where the radius rm has already disappeared inside the horizon, 
and is the result of the fact that at a fixed radius, decreasing 9 also decreases the 
flat space circular radius p = r sin 9, which increases the centripetal acceleration 
roughly by 1/p compared to the central attraction which docs not dramatically 
change with angle. This in turn increases the parallel transport rotation-dominated 
zone compared to the boost-dominated zone. The fact that the velocity profile 
corresponding to r 2 = (where I\ = 0), which determines the extremum of the 
acceleration polar angle \i is an extremely small negative function at these scales 
(explaining why all the I\ curves appear to pass through the origin in Fig. 2) shows 
that this point on the acceleration hyperbola is just slightly to the countcrrotating 
side of the ZAMO point. 

In the case of the equatorial plane of the Kerr spacetime discussed extensively 
in [15], the acceleration is purely radial corresponding to x = 0, 7r which makes 
its x-derivative and hence Ti identically zero. Thus A is always semi-singular and 
becomes singular at the boundary between boost and rotation dominance where 



Geometric transport along circular orbits in stationary axisyrametric spacetimes 27 



a) = 1.5 




Fig. 8. The figure corresponding to Fig. 6 off the equatorial plane (9 = 7r/2 fs 1.57) nearby at 
8 = 1.5 and farther away at 9 = n/3. Vanishing ti occurs at the solid curves (not bounding the 
grey region), while the dashed curve indicates the zero of its ^-derivative. In a) one sees that 
departing from the equatorial plane, the geo— and counterotating crit— curves to the right of rj 
in the equatorial plane Fig. 6 deform and join together to produce the subluminal counterrotating 
vanishing t\ curve, while the geo+ curve and the corotating crit— curve similarly produce the 
corotating such subluminal curve; the counterrotating superluminal crit+ curve joins together with 
the superluminal geo— curve to produce the the superluminal counterrotating vanishing n curve. 
The boost zone shrinks with decreasing angle. The dashed-dotted uppermost curve in the white 
rotation-dominated zone corresponds to vanishing T2 at positive superluminal speeds (converging 
with the boost-dominated zone boundary curves at the horizon at v = 0), while the subluminal 
such curve occurs at very small negative relative velocities not distinguishable from the i/-axis at 
this scale. 
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\k\ — \t\\. This occurs for 

"(gmp) T V+ V- 



/ {PT±) 



1 T "(gmp) 



f(gmp) = {v+ +V-)/2 , (73) 



where ^( gm p) 1S the velocity of the geodesic meeting point trajectories, 15 or explicitly 
to a{r-m)(a 2 + 3 r 2 ) ± (r 3 + a 2 r + 2 a 2 m) \fK 



V(PT±) 



(74) 



(— r 4 + 2 m r 3 + 2 a 2 to r + a 2 m 2 ) \/A 

Notice that from the formula (73), the values the values vr PT ±\ = ±1 occur when 
respectively v± = ±1, so that the respective endpoints of the boost-dominated 
zone go null at the same radii as the counter-rotating (rj) and then the corotating 
circular geodesies (ru) go null as one approaches the horizon. Similarly, for example, 
as ^( g mp) ~> ~ 1 + j then ^ PT _j — > —00, which occurs at rm in Fig. 6. This feature 
of the geodesic meeting point horizon in the equatorial plane was uncovered in 
studies of circular holonomy started by Rothman et al. 15-17 

From Eqs. (4.9) and (4.10) of [13], the conditions v± = ±1 are equivalent to 
^(crit)± = ±1; so m turn the conditions V(pt-) = — 1 an d ^(pt+) = 1 correspond to 
the horizons 77 and rjj of the extremely accelerated orbits in the equatorial plane, 
which are critical points of the acceleration as a function of the relative velocity, 
labeled as crit± in Fig. 6. This explains the triple crossing points of the geodesic 
velocities, extremely accelerated velocities and the boost zone endpoint velocities 
in that figure. The solutions of the conditions ^(pt±) = ±1, when fractions are 
cleared (numerator equals plus or minus the denominator) to produce polynomial 
equations in r after some manipulation to eliminate the radical, have a common 
factor 4 a 2 m — 9 m 2 r + 6 m r 2 — r 3 , whose zeros produce the radii rj and ru of the 
pair of null circular geodesies. The zero of the denominator of Eq. (74) is instead 
rm — f( g mp)7 the horizon of the geodesic meeting point observers. 

As a/m increases from the value 1/2 of Fig. 6 to the extreme limit a/m = 1 in 
the equatorial plane, the radius ru shrinks down to the horizon r( hor ) crossing over 
rm as shown in Fig. 2 of 15 where these radii are labeled respectively by "geo+," 
"gmp" and "hor." In fact very interesting changes occur in the interval of values 
[0.89, 1] for a/m as illustrated in the polar plots of Fig. 9. As the surface II shrinks 
down to the horizon, the surfaces 177 and I both bubble upwards and outwards, 
with surface / developing a transversal intersection with the horizon only at the 
limiting value a/m — 1 of an extreme Kerr spacetime. That surface i7 shrinks to the 
horizon in this limit means that there always exists a rotation-dominated zone for 
sufficiently fast but subluminal corotating orbits outside the horizon. The increasing 
separation between the horizon and surface III away from the North pole in this 
limit means that the corresponding small piece of grey boost-dominated region at 
the top left corner of Figs. 6 and 8a) grows outward and downward as the lower grey 
zone falls downward following the rotation-dominated region as its lower border 
in turn moves downwards. Simultaneously the counterrotating boost-dominated 
region moves outward from the horizon away from the North pole. Thus the faster 
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Fig. 9. The figure corresponding to Fig. 7 as a/m increases from 1/4 towards the extreme value 
1. The surface II collapses to the horizon, while the surfaces I and III bubble upwards and 
outwards (like the ergosphere boundary which they cross) to form two doughnut-likc rings (when 
revolved about the z-axis) around the horizon, one inside the other. In the Schwarzschild limit 
a/m — ► 0, the surfaces I and 77 collapse together to a single surface meeting the equatorial plane 
at r/m = 3, while the ergosphere boundary and surface III collapse together to the horizon. 

the spacetime rotates, the worse things are for the counterrotating orbits, which go 
superluminal farther and farther away from the shrinking horizon leading to boost- 
domination inside that radius for all subluminal motion. At the same time the 
corotating orbits go superluminal closer and closer to the horizon, with expanding 
rotating-dominated zones existing for suitably fast but subluminal motion outside 
that radius. 

8. Concluding remarks 

The gravitoclcctromagnetic analysis of parallel transport of the 4-velocity along 
generic circular orbits in stationary axisymmctric spacetimes has been extended 
to parallel transport of the full tangent space and connected to the corresponding 
Serret-Frenet description. The electromagnetic decomposition of the parallel trans- 
port matrix is essential to putting it into a canonical form for interpretation. Crucial 
features of parallel transport take place at the surfaces where various geometrically 
defined 4-velocities become null. These results enable the previous discussion of 
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the equatorial plane of the Kerr spacetime to be extended to the whole exterior 
spacetime outside the horizon. 

One finds new surfaces surrounding a black hole which intrinsically characterize 
limiting parallel transport properties along circular orbits of all causality types. 
Within these surfaces, one has an invariant characterization of the domination of 
centripetal acceleration effects by attractive forces, differing of course for corotating 
and counterrotating orbits. Examining the way in which the gravitoelectric, gravit- 
omagnetic and space curvature effects enter into the parallel transport matrix and 
the resulting 4-acceleration gives some more insight into how centripetal accelera- 
tion and gravitational attraction compete with each other and are in turn affected 
by gravitomagnetic effects in the Kerr spacetime, thus complementing the results 
of previous work. 9 



Appendix A. Geometry of the acceleration and Prenet-Serret 
angular velocity curves 

Using the defining condition for the MIROs (vertex observers), (see the paragraph 
containing equation (6.3) of [9]) e 2a < vort > = y/\ |a_ 1 1/| \a + \ |, and introducing the 
quantity A = ^/||a_||||a+||/2, the definitions of section 4 together with the decom- 
position a± = \ \a±\\n± into lengths and unit direction vectors lead to 

{3 -So) = \ [||a + ||e 2Q n + + ||a_|| e - 2Q n_] 

= ^ [ e 2 («-«(vert)) n + + e -2(a-a (vert) ) ( A ^ 

A 

= — [cosh2(a - a(vert)) (n+ + n_) + sinh2(a - a( ver t)) (n+ - ft-)] . 

Here a = arctanhj/ is the rapidity paramctrization of the acceleration curve. 

The photon acceleration direction vectors, which determine the common asymp- 
totes of the acceleration and Frenet-Serret angular velocity hyperbolas, are inde- 
pendent of the reference observer family used to decompose the acceleration. It is 
natural to introduce the unit direction vectors of their sum and difference which 
determine the common orthogonal axes of the acceleration hyperbola. 

The opening angle j3 £ [0, n] of the acceleration hyperbola (making n — (3 the 
opening angle of the Frenet-Serret angular velocity hyperbola) satisfies ri + ■ n_ = 
cos (3, in terms of which one finds 

||n+ + n_||/2 = cos^ , ||n+ - n_||/2 = sin | , (A.2) 
leading to the orthogonal unit vectors 

r= £±£ • r- (A.3) 

2 cos 77 2 sin ^ 

which give the directions of the principal axes of both hyperbolas. The quantities 
A cos § and A sin | are the semi-axis parameter lengths of the acceleration and the 



Geometric transport along circular orbits in stationary axisymmetric spacetimes 31 



Frenet-Serret angular velocity hyperbolas. By introducing the sign e which makes 
the following cross-product relation valid 

eix n j = e f x„ e § = , (A.4) 

and whose value in the Kerr spacetime is e = sgn(7r/2 — 6), one can finally write 
the equations for the acceleration and angular velocity hyperbolas in terms of their 
geometric properties 

a - a = .4[cosh2(a - a (vert )) cos ^ f+ sinh2(a - a (vert )) sin ^ j] , 

e^(FS) = ^4[cosh2(a - a (vcrt )) sin ^ T- sinh2(a - a (vcrt) ) cos^j] . (A.5) 

It is then easy to derive the following formulas for the acceleration and angular 
velocity scalars 

2 (a - a ) • W( FS ) = A 2 sin (3 , 

\\a — a || 2 = ^ ^4 2 [cosh4(a — "(vert)) + cos/3] , 
ll^(FS)H 2 = - ^ 2 [cosh 4(a - a (vcrt) ) - cos /3] , (A.6) 

or equivalently 

\\a — a \\ 2 - ||cj( FS )|| 2 = A 2 cos/3 , 

||a-a || 2 + ||u; (FS) || 2 = A 2 cosh4 (a - a (vcrt) ) . (A.7) 



Appendix B. Frenet-Serret frame for circular orbits 

The sign conventions for the Frenet-Serret scalars and the corresponding directional 
choices of the frame vectors for nonnull circular orbits in stationary axisymmetric 
spacetimes require careful consideration. 9 Starting from the rapidity parametrized 
4- velocity for a future-pointing timelike circular orbit 

U = coshan + sinhae^ = e (B.l) 

expressed in a Boyer-Lindquist-like ZAMO orthonormal frame {n, ef, eg, e^}, one 
can introduce a unit vector in the direction of increasing </> in the local rest space 
of U, namely 

d 

e 2 = — — U = sinh a n + cosh a e ? , (B.2) 

da v 

thus defining two Frenet-Serret vectors in the circular orbit velocity plane so that 
eo A e 2 = n A e It remains to be seen that e 2 is actually a possible choice. 

Next one introduces polar coordinates in the acceleration plane spanned by 
{e r , eg} to parametrize the acceleration vector, with a nonnegative radial coordinate 

a(U) = k (coax ef + sinx eg) — nei , K >0, (B.3) 
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so that ei is the unit direction of a(U) and e 3 is the orthogonal unit vector 

d . . 

e 3 = -— ei = sin%er - cos%e e - , (B.4) 

satisfying 

e^x„ei = -e 3 , x„ e 3 = ei , e 3 A ei = A , (B.5) 

and so eo A ei A e2 A e 3 = n A A eg A e^. Since n and eo both future-pointing 
timclike vectors, the Frenet-Serret spatial triad is also right-handed. 

To show that e 2 and e 3 can consistently be chosen as the remaining Frenet- 
Serret vectors, the transformation (39) of the connection matrix A must be used 
together with the identifications (42), leading to 

T 1 )= M ^( g n) + ( 9 a 

T 2 e 1 +T 1 e 3 J \H ) V 



(B.6) 



(cosh 2a g- + sinh 2a x „ H + g + 
— sinh 2a x „ + cosh 2a H 

Then since 5+ is independent of a and x„ (e^ x„ i7) = —H, the Frenet-Serret 
angular velocity can be rewritten as 

1 d 1 dn 1 d\ dei 

r 2 e 1 +r ie3 = - e ^x„-(.e 1 )=2^e^x„e 1 + -.- e ^x rl — 

= -2^-2 K d^ ei - (R7) 
Thus everything is consistent if one makes the identification 

1 dn 1 dx ,„ e s 

ri = -2^' ^r^- (R8) 

The Frenet-Serret angular velocity is then given by the rapidity derivative formula 
W(FS) = ^ e x„ = cosh 2a # - sinh 2a x n g_ . (B.9) 

Note that interchanging the hyperbolic functions in (B.l) leads to the corotating 
spacelike case for U, which only interchanges U and resulting in the same iden- 
tifications (B.8) and the same first equation of (B.9). It remains to check all these 
signs for the remaining two cases of the four sign combinations for parametrizing 
U by a hyperbolic angle a as in Eq. (17). 



References 

1. B.R. Iyer and C.V. Vishveshwara, Phys. Rev. D48, 5706 (1993). 

2. R.D Greene, E.L. Schiicking and C.V. Vishveshwara, J. Math. Phys. 16, 153 (1975). 

3. C.W. Misner, K.S. Thorne and J. A. Wheeler, Gravitation, (San Francisco: Freeman), 
1973. 

4. F. de Felice, Class. Quantum Grav. 12, 1119 (1991); F. de Felice and S. Usseglio- 
Tomasset, Class. Quantum Grav. 8, 1871 (1991). 

5. D. Page, Class. Quant. Grav. 15, 1669 (1998). 



Geometric transport along circular orbits in stationary axisyrametric spacetimes 33 



6. O. Semerak, Gen. Relativ. Grav. 30, 1203 (1998). 

7. D. Bini, F. dc Felice and R.T. Jantzen, Centripetal acceleration and centrifugal force 
in general relativity in Nonlinear Gravitodynamics. The Lense-Thirring effect, Ed. R. 
Ruffini and C. Sigismondi (Singapore: World Scientific), 2003. 

8. D. Bini, F. de Felice and R.T. Jantzen, Class. Quantum Grav. 16, 2105 (1999). 

9. D. Bini, R.T. Jantzen and A. Merloni, Class. Quantum Grav. 16, 1333 (1999). 

10. J.L. Synge, Relativity: The General Theory, North Holland Amsterdam, 1960. 

11. R.T. Jantzen, P. Carini and D. Bini, Ann. Phys. (N.Y.), 215, 1 (1992). 

12. D. Bini, P. Carini and R.T. Jantzen, Int. J. Mod. Phys. D6, 1 (1997). 

13. D. Bini, P. Carini and R.T. Jantzen, Int. J. Mod. Phys. D6, 143 (1997). 

14. D. Bini, P. Carini and R.T. Jantzen, Class. Quantum Grav. 12, 2549 (1995). 

15. D. Bini, R.T. Jantzen and B. Mashhoon, Class. Quantum Grav. 19, 17 (2002). 

16. T. Rothman, G.F.R. Ellis and J. Murugan, Class. Quantum Grav. 18, 1217 (2001). 

17. R. Maartens, B. Mashhoon and D.R. Matravers, Class. Quantum Grav. 19, 195 
(2002). 

18. D. Bini, R.T. Jantzen and B. Mashhoon, Class. Quantum Grav. 18, 653 (2001). 

19. D. Bini, C. Cherubini and R.T. Jantzen, Class. Quantum Grav. 19, 5481 (2002). 

20. J.L. Synge, Relativity: The Special Theory (2nd Edition), Chapter IV, North Holland 
Amsterdam, 1965. 

21. C.G. Bollini, J.J. Giambiagi and J. Tiomno, Lett. Nuovo Cim. 31, 1 (1981). 

22. D. Bini, F. de Felice and R.T. Jantzen, Class. Quantum Grav. 16, 2105 (1999). 

23. D. Bini and R.T. Jantzen, Proceedings of the Ninth ICRA Network Workshop on Fermi 
and Astrophysics, eds. V. Gurzadyan and R. Rufiini, World Scientific, Singapore, 2003. 



